Application of a symmetry adapted algebraic 
model to the vibrational spectrum of methane 
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^ ! 1 Algebraic model 

Recently, a symmetry-adapted version of the vibron model has been developed [1- 
4]. The aim of this work is to present the description of all stretching and bending 
J>"v vibrations of the methane molecule in the framework of such a model. 

In the present algebraic approach a U(2) algebra is associated with each relevant 
interatomic interaction. For the CH4 molecule we have four U(2) algebras correspond- 
1 I ing to the C-H interactions and six more representing the H-H couplings. The assign- 

ments and the choice of the Cartesian coordinate system are the same as in []. The 
^ ■ molecular dynamical group is then given by the product Z7i(2) ® . . . (g) Ui (2), where 

each Ui(2) algebra (i — 1, . . . , 10) is generated by the set {Gi} = {N^ J+j, J- t i, Jo,i}, 
which satisfies the commutation relations 

: 

^ ' &i, = ±J±,i , [J+,i, J-A = 2J ,i , [Ni, 4i] = , (1) 

with fi = ±, 0. The labeling is such that i = 1, . . . , 4 correspond to the C-H couplings 
while the other values of i are associated with H-H interactions []. Here Ni is the 
i-th boson number operator. All physical operators are expressed in terms of the 
generators {Gi}, and hence commute with the number operators N{. Since Jf = 
.£h I Ni(Ni + 2)/4 we can make the identification ji = N/2. The eigenvalues of Jo.j are 

restricted to m, > and can have the values = N/2, (N — 2)/2, ...,1/2 or 
for Ni odd or even, respectively []. The local basis states for each oscillator are 
usually written as \N,Vi), where Vi = (N — 2rrii)/2 = 0, 1, . . . [N/2] denotes the 
number of oscillator quanta in the i-ih oscillator. For the CH 4 molecule there are 
two different boson numbers, N s for the C-H couplings and Nf, for the H-H couplings, 
which correspond to the stretching and bending modes, respectively. 

The symmetry of methane is taken into account by projecting the local operators 
{Gi}, which act on bond i, on the irreducible representations T of the tetrahedral 
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group %. For the J^i generators of Eq. (|I|) we obtain the Td tensors 

10 

f ln = E«S4^ (2) 

i=l 

where 7 denotes the component of T, and the label x refers to stretching (s) or 
bending (b). The explicit expressions for the tensors are similar to the one phonon 
wave functions expansions given in []. The algebraic Hamiltonian is constructed by 
repeated couplings of these tensors to a total symmetry A\. 

The methane molecule has nine vibrational degrees of freedom. Four of them 
correspond to the fundamental stretching modes (Ai © F 2 ) and the other five to the 
fundamental bending modes (E(BF 2 ) []. The projected tensors of Eq. (2) correspond 
to ten degrees of freedom, four of which (Ai © F 2 ) are related to stretching modes 
and six (A\ Q)EQ)F 2 ) to the bendings. Consequently we can identify the tensor T^\ b 
as the operator associated to a spurious mode. This identification makes it possible 
to eliminate the spurious states exactly. This can be achieved by (i) constructing the 
physical Hamiltonian by simply ignoring the T^l' b tensor in the coupling procedure, 
and (ii) diagonalizing this Hamiltonian in a symmetry-adapted basis from which the 
spurious mode has been removed following the procedure of []. 

According to the above procedure, we construct the 7^ invariant interactions that 
are at most quadratic in the generators and conserve the total number of quanta 

wr. = ^ E rfe + t% f%) , 

^ = ^E^ r ;- (3) 

Here T = A\, F 2 for the stretching vibrations x = s and T = E, F 2 for the bending 
vibrations x = b. In addition to Eq. (^), there are two stretching-bending interactions 
which will not be considered at the moment. 

The zeroth order vibrational Hamiltonian is then written as 

Hq = LOi Hai >s + ^2 T~(-E b + ^3 7~(-F2,s + ^4 'HF 2>b 

+a 2 V Eb + «3 Vf 2 , s + «4 Vf 2 , 6 • (4) 

The interaction Va 1s has not been included since the combination J2ri^r a + Vr s ) 
corresponds to the constant N s + 2. A similar situation is present for the bending 
interactions, but in this case the interaction b has already been excluded in order 
to remove the spurious A\ bending mode. The subscripts of the parameters corre- 
spond to the (y\, v l 2 , z/3 3 , z/4 4 ) labeling of a set of basis states for the vibrational levels 
of CH 4 . Here u%, u 2 , ^3 and denote the number of quanta in the A\ tS , Eb, F 2yS and 
F 2t b modes, respectively. The labels U are related to the vibrational angular momen- 
tum associated with degenerate vibrations []. In the harmonic limit the interactions 



of Eqs. (|3]) and (4) attain a particularly simple form, which can be directly related 
to configuration space interactions [1-4]. The / Hr x terms represent the anharmonic 
counterpart of the harmonic interactions, while the Vr^ terms are purely anharmonic 
contributions which vanish in the harmonic limit. 

The zeroth order Hamiltonian of Eq. (|J) is not sufficient to obtain a high-quality 
fit of the vibrations of methane. Several other physically meaningful interaction are 
essential for such a fit. For the study of the vibrational excitations of methane we 
use the 7^ invariant Hamiltonian 

H — ui H Al s + u 2 H Eb + u 3 H F2 a + u A H F2 b + a 3 V F2 s 

-hXn (j~tAi, s ) + X 2 2 (ftE b ) + X 33 (Wf 2 , 8 ) + X u {Hf 2 ^ 
+X12 (fCA hs HE b ) + Xu \HAi,s ^F 2ib ) 

+X23 (Ti.E b / Hf 2 , s ) + X24 {jrf-E h ^F 2i6 ) + -^34 (i~iF 2 , s ^#2,&) 

+922 + 933 E ^ ?7 + 9* £ ^ ^ + 934 E 

7 7 7 

+^33 @ss + ^44 Cfefe + ^34 Osb + ^23 @2s + ^24 2 b ■ (5) 

The Xij terms are quadratic in the operators 1tir x and hence represent anharmonic 
vibrational interactions. The terms are related to the vibrational angular momenta 
associated with the degenerate vibrations [1-3]. 

In the harmonic limit the expectation value of the diagonal terms in Eq. (^) leads 
to the familiar Dunham expansion []. 

I>(«i+*)+EE (*« + % + y ) + E E 9a kh • (6) 

i j>i i j>i i 

Here di is the degeneracy of the vibration. The £y terms in Eq. (|5|) give rise to further 
splittings of the vibrational levels (z/ 1; z/ 2 , z/ 3 , z/ 4 ) into its possible sublevels []. Their 
explicit expressions are given in ref. [2]. 

2 Results 

The Hamiltonian of Eq. (|5]) involves 23 interaction strengths and the two boson 
numbers, N s and iV&. The vibron number associated with the stretching vibrations 
is determined from the spectroscopic constants u e and x e u e for the CH molecule to 
be N s = 43 []. The vibron number for the bending vibrations, which are far more 
harmonic than the stretching vibrations, is taken to be Nf, = 150. We have carried 
out a least-square fit to the vibrational spectrum of methane including 44 energies. 
We find an overall fit to the observed levels with a r.m.s. deviation which is an order 
of magnitude better than in previous studies. While the r.m.s. deviations of [] and [] 
are 12.16 and 11.61 cm -1 for 19 energies, we find a r.m.s. of 1.16 cm -1 for 44 energies. 



The values of the fitted parameters as well as all predicted levels up to V = 3 are 
given in ref. [1]. 

The «3 term is completely anharmonic in origin and has no counterpart in the 
harmonic limit. In order to address the importance of this term in Eq. ([|) we have 
carried out another calculation without this term. With one less interaction term 
the r.m.s. deviation increases from 1.16 to 4.48 cm -1 . This shows the importance of 
the term proportional to 0:3 to obtain an accurate description of the anharmonicities 
that are present in the data. The absence of the a.3 term in the second calculation 
can only partially be compensated by the anharmonicity constants Xij []. 

3 Conclusions 

In summary, in this work we present the description of vibrational excitations of 
methane in a symmetry-adapted algebraic model. We find an overall fit to the 44 
observed levels with a r.m.s. deviation of 1.16 cm -1 , which can be considered of 
spectroscopic quality. We pointed out that the Vf 2s term in combination with the 
anharmonic effects in the other interation terms plays a crucial role in obtaining a fit 
of this quality. Purely anharmonic terms of this sort arise naturally in the symmetry- 
adapted algebraic model, but vanish in the harmonic limit. Physically, these con- 
tributions arise from the anharmonic character of the interatomic interactions, and 
seem to play an important role when dealing with molecular anharmonicities. 

These studies suggest that the symmetry-adapted algebraic model provides a nu- 
merically efficient tool to study molecular vibrations with high precision. The main 
difference with other algebraic methods is the use of symmetry-adapted tensors in 
the construction of the Hamiltonian. In this approach, the interactions can be con- 
structed in a systematic way, each term has a direct physical interpretation, and 
spurious modes can be eliminated exactly. 
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